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The structural properties of fluids whose molecules interact via potentials with a hard core plus two piece- wise 
constant sections of different widths and heights are presented. These follow from the more general devel- 
opment previously introduced for potentials with a hard core plus n piece-wise constant sections [Condons. 
Matter Phys. 15, 23602 (2012)] in which use was made of a semi-analytic rational-function approximation 
method. The results of illustrative cases comprising eight different combinations of wells and shoulders are 
compared both with simulation data and with those that follow from the numerical solution of the Percus- 
Ycvick and hypernettcd-chain integral equations. It is found that the rational-function approximation predicts 
a more accurate radial distribution function than the Percus-Yevick theory in all the cases and is comparable 
to or superior than the hypernetted-chain theory in most of the cases. 



I. INTRODUCTION 

Due to their relative simplicity, while having at the 
same time the ability to adequately account for diverse 
physical features of real fluids, discrete potentials of the 
form 



ip{r) 



oo 



r < a, 

(J < r < Aicr, 



£2, Aio- < r < X^cr, 



e„, A„_ifT < r < A„CT, 
, 0, r > A„o-, 



(1) 



have received some attention in the recent literatureii"— 
They comprise a hard core of diameter a and n steps of 
"heights" Cj and widths {Xj — Aj_i )cr, with Aq = 1, so 
that A„cr denotes the total range of (^(r). The sign of ej 
determines whether the jth step is either a "shoulder" 
{ej > 0) or a "well" {ej < 0). The interaction potential 
at r = Xjcr {j = l,2,...,n) is repulsive if ej > ej_|-i 
and attractive if ej < ej^i (with the convention e„-|-i = 
0). Particular cases of these discrete potentials when 
n — \ are the popular square-well and square-shoulder 
potentials. 



^'Electronic mail: 

|http://www.unex.es/eweb/fisteor/andres7] 
"•' Electronic mail: 

mttp : / /ww w. unex.es/eweb/fisteor/santos/ 1 

'^' Electronic mail: 

[http://xml.cie.unam.mx/xml/t c/ft/mlh/| 
'''Electronic mail: mbarccnas@ tese.edu. mxl 
*' Electronic mail: jporea@imp.nix] 



landres@unex.esl 

lsantos@unex.esl 

|malopez@unam.mx[ 



The phase diagram and the thermodynamic proper- 
ties of discrete-potential fluids have been thoroughly ex- 
amined and are relatively well understood4r"— On the 
other hand, studies of their structural properties either 
theoretical or from simulations are scarceJ^ii^r— In a 
previous paper fi^ following a semi-analytic methodology 
referred to generically as the rational-function approxi- 
mation (RFA) that, although approximate, has proven 
successful for many other systems?^ we derived the gen- 
eral formulae for the structural properties of fluids whose 
molecules interact via discrete potentials with a hard core 
plus an arbitrary number of piece- wise constant sections 
of different widths and heights. The theoretical scheme 
was illustrated by comparing it with available computer 
simulations resultsj^^ 

The aim of this paper is to carry out a more systematic 
study of the structural properties of fluids characterized 
by a discrete potential with a hard core plus different 
combinations of a repulsive shoulder and an attractive 
well. This will be done by considering the results of Rcf. 
Il9l in the case of n = 2 for various values of the param- 
eters and subsequently performing a comparison both 
with simulation results as well as with those that fol- 
low from the numerical solution of the Ornstein-Zernike 
(OZ) equation with both the Percus-Yevick (PY) and 
hypernetted-chain (HNC) closures. As will be seen, the 
performance of the RFA approach, despite its simplicity, 
is quite satisfactory. 

The paper is organized as follows. In order to make it 
self-contained, in Sec. lUwe introduce the systems to be 
studied and sketch the derivation of the results of Ref. [l^ 
for the structural properties of such systems when n = 2. 
We also include here some details of the simulation and 
of the numerical solution of the OZ equation with the 
PY and HNC closures. This is followed in SecHlIlby the 



comparison between the outcomes of the three different 
approaches for the radial distribution function (rdf ) . The 
paper is closed in Sec. IIVI with further discussion and 
some concluding remarks. 



II. SYSTEM AND STRUCTURAL PROPERTIES 

We consider a fluid of number density p and absolute 
temperature T in which the intermolecular pair potential 
is of the form of Eq. ([T|) with n = 2. We will take the 
hard-core diameter a as the length unit so all distances 
will be measured in units of a. The two main quantities 
usually employed to characterize the structure of fluids in 
equilibrium are the static structure factor S{q) and the 
rdf g{r) which arc related by 



S{q) = l + p / dre-'^^-'lgir)-!] 
G{s)-G[-s) 



1 -27rp 



s=tq 



where 



Gis) 



dr e ^^rg{r) 



(2) 



(3) 



is the Laplace transform of rg{r). 



A. The rational-function approximation method 

We define an auxiliary function F[s) directly related 
to G{s) through 



G{s) 



F{s)e- 



1 + 12r;F(s)e-'' 



(4) 



Here, rj ~ {Tr/6)pa^ is the packing fraction. Laplace in- 
version of Eq. (d]) provides a useful representation of ^(r), 
namely 

oo 

g{r) = r-i ^ (-12r;)'"- V„,(r - m)e{r - m), (5) 

m— 1 

where fm{r) is the inverse Laplace transform of s[F{s)]"^ 
and 8(r) is the Heavisidc step function. 

The contact value g(l^) of the rdf is related to F{s) 
through ,9(1^) = ./i(0) = hm s~y 00 s'^F{s) and it has to be 
finite. Further, as seen from Eq. ([2]), the behavior of G{s) 
for small s determines the value of 5(0), which must also 
be finite. Hence, F{s) must satisfy two conditionsJ^ 



Fis) 






S — > CXD, 



l + 2r/ 



+ V 4 



12?/ 



24r; 



(6) 



(7) 



Equations (|4|)-([2]) are exact and valid for any inter- 
action potential with a hard core at r = cr = 1. Now 
we particularize to the potential ([T]) with n = 2. To re- 
fiect the discontinuities of g(r) at the points r = Ai and 
r = A2, where ip(r) is discontinuous, we decompose F{s) 
as 

F{s) = Ro{s) + i?i(s)e-(^i-i)-' + i^2(s)e-(^^-l)■^ (8) 

As a consequence, 

/i W = UrMr) + Ur - Ai + l)e(r - Ai + 1) 



+6(^-A2 + l)e(r-A2 + l), 



(9) 



where S.j{r) denotes the inverse Laplace transform of 
sRj{s). If, as will be done here, one assumes that A2 < 2, 
insertion of Eq. ([9]) into Eq. ([5]) gives the rdf in the shell 
1 < r < 2. In particular. 



ff(Ar) = Ar^eo(Ai - 1), 



5(A^) = ArMeo(Ai-l)+6( 



(10) 



(11) 



5(A2-) = A^-i [eo(A2 - 1) + ei(A2 - Ai)] , (12) 

ff(Aj) = A^-i Ko(A2 - 1) + 6(A2 - Ai) + 6(0)] . (13) 

Now we assume the following rational-function approx- 
imation for Rj{s): 

R.(s)--^ ^^+^^' 7-012 

(14) 
Note that Eq. p^ for Ro{s) guarantees the fulfillment of 
the physical condition (|6]). Yet, the approximation p4)) 
contains nine parameters to be determined. The exact 
expansion ^ imposes five constraints among those nine 
parameters, nameljii^ 



1 = Ao + Ai+A 



Si 



-l + Bo-C^'\ 



S3 



S2 = l-Bo + C^'^ + lc^'\ 



1 + 2'? , 1,_ 1^(1) 1^(2) 1 



(15) 



(16) 



(17) 



12, +^^°-^^'^'-^^'^'-^^'"' (^«) 



B, = C(i) + ^^ f6C(2) + 4C(3) + CW) + Hull, 
1 + 2?/ V / 1 + 277 



(19) 



Here, 






Next, since the cavity function y{r) = g{r)e^'^^'^\ where 
/3 = l/k^T (kg being the Bohzmann constant), must be 
continuous at r = Ai and ?' = A2, one obtains from Eqs. 
(Unil-dlSl) the two conditionsiS 



S3 



El 

S3 



,/3(ei-e2) 



'IE 



,(Ai-l)s„ 



^ 5l + 2S2S0, + 3S3SI 

Q — 1 " 



x{Ao+BoSa), 



(21) 



„/3e2 



An 



3 



Srx€ 



(A2-l)s„ 



^ Si + 2S2Sa + 35*352 

a— 1 



,(22) 



where Sq, (a ~ 1,2,3) are the three roots of the cubic 
equation 

1 + SiS^ + S2SI + ^34 = 0. 



(23) 



Equations ([151) -(O, (ED, and ^ still leave two pa- 
rameters undetermined. A simplifying assumption is that 
the coefficients Aj [j = 0, 1, 2) may be fixed at their zero- 
density values, namely 



^0 



-/Jei 



Al 



-/3ei 



1 



-/3e2 



(24) 

This closes the problem of determining the nine param- 
eters in terms of 77, Ai, A2, /3ei, and j5e2- In fact, Eqs. 
PB)) ^ ([TO)) allow us to express Si, S2, S3, and Bq as linear 
combinations of Bi and B2, so that in the end one only 
has to solve (numerically) the two coupled transcendental 
equations (PT|) and (P^ . Since the dependence of G{s) 
on s is explicit, we are now in a position to compute the 
structural quantities of our systems. The structure fac- 
tor S{q) can be directly obtained from Eq. ([2]), while the 
rdf g{r) can be obtained from Eq. ([5|) or, more directly, 
by numerical inverse Laplace transform of G{s)i^ 



B. The PY and HNC approximations 

In the usual integral equation approach to the theory 
of liquids, the OZ equation, which may be formally con- 
sidered as a definition of the direct correlation function 
c{r) , provides a link between this direct correlation func- 
tion and the total correlation function h{r) = g{r) — 1, 
the latter being a measure of the 'influence', either direct 
or through a third molecule, of two molecules separated 
a distance r away. The OZ relation reads 



h{r) = c{r) +p dr'c{r')h{\r - r'|) 



i-) />oo pr-\-r 

^ c{r) -] / dr'r'c{r') I dr"r"h{r"), 

1" Jo J\r-r'\ 

(25) 

where in the second equality we have particular- 
ized to three-dimensional systems and used bipolar 
coordinates i^i^ 



Since both h{r) and c{r) are unknown, in order to 
close the description one requires an additional equa- 
tion, known as the closure relation. A closure can be 
expressed as a local relationship between the direct corre- 
lation function and the Mayer function /(r) = e"^''''-''-' — 1 
and the cavity function y{r), i.e.. 



c{r)=Cifir),y{r)) 



(26) 



Equivalently, Eq. (pS)) can be inverted to obtain a local 
relationship between the cavity function and the indirect 
correlation function j{r) = h{r) — c{r), i.e.. 



y{r)^y{-l{r)). 



(27) 



Insertion of the closure ((26| and \27\ into the OZ relation 
(j25p yields a closed nonlinear integral equation for the 
cavity function: 



y{r) = y 



-(^. 

r-\-r' 



X / dr r' 

I |j-_^' I 



dT'r'C{f{r'),y{r')) 
-0^(^-")y{r") - 



(28) 



As mentioned in Sec. HI we will consider here both the 
PY and HNC closures given by 



c{r) = h{r) - y{r) + 1 (PY), 



c{r) ^h{r)~lnyir) (HNC). 



(29) 



(30) 



In terms of the functions C{f,y) and 3^(7), the PY and 
HNC closure relations are 



Cif,y)^fy, yii) = i + i (PY), 



(31) 



C{f,y) = {f+l)y-l-\ny, 3^(7) = e^ (HNC). (32) 



We have solved Eq. ([^5)) numerically in the PY and 
HNC cases. First, a discretization scheme y{r) — > {r^, y,} 
with Ti = iAr (i = 1,2,..., A/") and a cut-off distance 
r_\f = A/'Ar are introduced, so that the integral equation 
(PS]) is replaced by a set of M nonlinear coupled equa- 
tions: 



y^ = y 



2ttp 



Af 



(Ar)2^r,C(/„y,) 



j=i 



niin{i-\-j.jV) 

X Y. ^fc (e-^'^'y^ - 1) 1 , ^^l,2,...,^f. 

k=\i~j\ + l 

(33) 

Next, a coarse-grained solution of Eq. ((55)) is obtained 
by an iteration method. Insertion of the nth order input 
{2/1 "} into the right-hand side of Eq. (p3l) gives the 
nth order output {y, } and the subsequent input is 

; ; J (n+l,in) (n.out) /-, n (n,in) .,, 

constructed as y^ = ayl -I- (1 — ajyl , with 



a convenient choice of the mixing parameter a. The it- 
erations are continued until the convergence criterion 



I (n.out) (n.in) i ^ -, ri 

max]?/) -y, 'I < 10 



-3 



(34) 



is reached or the number of iterations exceeds 100. Once 
the coarse-grained solution has been obtained, a fine- 
grained solution of the set of Af equations (|33|) is de- 
termined with the help of a computational software pro- 
gram, using the coarse-grained solution as a seed<^ In the 
numerical solutions presented in Sec. IIIII we have gener- 
ally used Ar = O.Olcr and JV = 400. As for the mixing 
parameter a, it was chosen by trial and error, being nec- 
essary to decrease it as the density increases. 



D 



C B 



-H \ 1 



H 



C. Technical simulation details 

The simulation data were computed with a Replica 
Exchange Monte Carlo (REMC) method. The REMC 
method is also known as Parallel Tempering and was de- 
rived to achieve good sampling of systems that present 
a free energy landscape with many local minimal^ The 
REMC method consists of simulating Af replicas (copies) 
of the system at different thermodynamic conditions; the 
attempted MC moves are accepted or rejected according 
to the traditional Metropolis algorithm. Due to these ex- 
changes, a particular replica travels through many tem- 
peratures, allowing it to overcome any barriers to free 
energy. 

The method samples an expanded canonical ensemble, 
taking the temperature as the expansion variable. The 
existence of this expanded ensemble justifies the intro- 
duction of movements of exchange between replicas. The 
expanded ensemble is defined as 



Q. 



expanded 



M 



(35) 



where Qi is the partition function of the (NVT) canonical 
ensemble of the system (subensemble i) at temperature 
T, volume V, and number of particles iV. To satisfy the 
detailed balance condition, the probability of acceptance 
of the exchange is given by 



Pace = min(l,cxp[-(/3, - /3,)(t/, - U,)]), 



(36) 



where /3j — Pi is the difference between the reciprocal 
temperatures and Uj — Ui is the difference between the 
potential energies of the subensemblcs i and j. 

A cubic simulation box of dimensions L^ ~ Ly ^ L, ^ 
10(7 was used and periodic boundary conditions were set 
in the three directions. Verlct lists^Si were implemented 
to improve performance. We have carried out computer 
experiments for different systems, corresponding to dif- 
ferent values of the parameters of the potentials that will 
be specified later. These systems are further character- 
ized by their reduced density p* = pa^ and their reduced 



FIG. 1. (Color online). Diagram in the e2~e\ plane illus- 
trating the choices of the values corresponding to the eight 
different cases to be studied, labeled from A to H. 



temperature T* ~ ksTje. where e = max(|ei|, |e2|). 
The number of replicas M = 12 was chosen to match 
the number of different temperatures in which we want 
to examine the systems. The highest temperature was 
set at r* = 2, while the other temperatures were es- 
tablished following a decreasing geometric progression, 
namely T,t = 2(0.959)"-\ n = 1,...,12. The initial 
configuration of each system, consisting of a collection of 
A^ ~ 500 particles randomly arranged in the simulation 
box and thus setting the reduced number density of our 
systems as p* = 0.5, was equilibrated by conducting 10^ 
Monte Carlo simulation steps. The rdf was calculated 
over additional 4 x 10^ configurations. 



III. RESULTS 

For convenience and in order to try to be systematic, 
we now first fix the values of Ai and A2 to be Ai = 1.25 
and A2 = 1.5. As for the values of ei and £2 we have 
considered the eight representative cases sketched in Fig. 
[1] At this stage some comments are in order. Since, 
as already stated, e = max(|ei|, |e2|), at least one of the 
\ei\ must be equal to e. The other energy level has been 
chosen as 0, ±e, or zt^e with the following conditions. 
First, the cases having ei = £2 have not been considered 
since they correspond to having just one step. The same 
can be said about the cases with e2 = 0. Finally, if ei 
and £2 have opposite signs, then we have taken |ei| == 
|e2| = (-. The potentials of the different cases (from A to 
H) are represented graphically in Fig.[2l We observe that 
system A corresponds to a purely repulsive potential. In 
the cases B-D the potential is repulsive at r/cr ~ A2 
and attractive at r/a = Ai, with increasing attraction 
when going from B to D. System E presents a purely 
attractive tail beyond the hard core. As for systems F- 
H, the potential is attractive at r/a = A2 and repulsive at 
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FIG. 2. (Color online). Diagrams of the surveyed potentials, 
labeled from A to H. 



r/a — Ai, with increasing repulsion when going from F 
to H. The cases A and F-H are examples of core-softened 
potentials. 

As an illustration of the results of our calculations, 
and in accordance with the simulation experiments men- 
tioned above, we fix the reduced density to be p* = 
0.5 and consider the lowest reduced temperature T* = 
1.26193 (corresponding to n = M = 12) in all cases, al- 
though, as said before, we obtained simulation data for 
the whole temperature range 1.26193 < T* < 2. The re- 
sults for the rdf are displayed in Figs.lMTOlfor cases A-H, 
respectively. In order to assess the influence of ei and 62 
on g{r) at fixed (reduced) temperature and density, com- 
mon horizontal and vertical scales have been chosen in 
Figs. I3HT0I In the case of the purely repulsive system 
A, there exists a local accumulation of particles at the 
external edge (A^) of each repulsive step, followed by a 
local depiction at the internal edge (A~). For systems 
B-D, where the potential is repulsive at r/a = A2 but 
attractive at r/a = Ai, the population of particles (as 
seen by a reference particle at the origin) is depleted in 
the region Ai < r/a < A2 and increases when going from 
A^ to Aj~, as expected. These effects are enhanced as the 
depth of the inner well increases. In the case of system E 
the potential outside the hard core is purely attractive, 
what is reflected in an increase of g{r) from the exter- 



bo 




FIG. 3. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system corresponding 
to case A (Ai = 1.25, A2 — 1.5, ei = e, €2 — e/2) at p* — 0.5 
and T* = 1.26193 with simulation results. Solid line: RFA; 
dashed line: HNC; dotted-dashed line: PY; symbols: REMC 
data. The inset shows the interaction potential and the tem- 
perature value. 
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FIG. 4. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system corresponding 
to case B (Ai = 1.25, A2 = 1.5, ei = e/2, £2 = e) at p* = 0.5 
and T* = 1.26193 with simulation results. Solid line: RFA; 
dashed line: HNC; dotted-dashed line: PY; symbols: REMC 
data. The inset shows the interaction potential and the tem- 
perature value. 



nal edge A^ to the internal edge A^. Finally, systems 
F-H are the counterparts of systems B-D. Now the at- 
traction at r/a ~ A2 produces an increase of particles in 
the region Ai < r/a < A2, this effect being enhanced as 
the inner barrier becomes more repulsive. Interestingly 
enough, the correlation length, i.e., the distance beyond 
which \g{r) — 1| is lower than a certain small value, in- 
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FIG. 5. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system corresponding 
to case C (Ai = 1.25, A2 = 1.5, ei = 0, £2 = e) at p* — 0.5 
and T* — 1.26193 with simulation results. Solid line: RFA; 
dashed Une: HNC; dotted-dashed line: PY; symbols: REMC 
data. The inset shows the interaction potential and the tem- 
perature value. 



FIG. 7. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system correspond- 
ing to case E (Ai — 1.25, A2 — 1.5, ei = —e, 62 — — e/2) 
at p* — 0.5 and T* = 1.26193 with simulation results. Solid 
line: RFA; dashed line: HNC; dotted-dashed line: PY; sym- 
bols: REMC data. The inset shows the interaction potential 
and the temperature value. 
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FIG. 6. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system corresponding 
to case D (Ai — 1.25, A2 = 1.5, ei — — e, 62 = e) at p* = 0.5 
and T* — 1.26193 with simulation results. Solid line: RFA; 
dashed Une: HNC; dotted-dashed line: PY; symbols: REMC 
data. The inset shows the interaction potential and the tem- 
perature value. 



FIG. 8. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system correspond- 
ing to case F (Ai — 1.25, A2 = 1.5, ei = — e/2, £2 = — e) 
at p* = 0.5 and T* = 1.26193 with simulation results. Solid 
line: RFA; dashed line: HNC; dotted-dashed line: PY; sym- 
bols: REMC data. The inset shows the interaction potential 
and the temperature value. 



creases monotonically from system A (Fig. [3]) to system 
H(Fig.[lOl). 

Let us comment now on the theoretical predictions. It 
follows that, as already pointed out in Ref. [l^ the RFA 
approach certainly outperforms the PY approximation 
in all the cases. This is noteworthy because, while the 
RFA reduces to the PY solution for hard spheres^i^ it is 



much simpler than the PY integral equation theory for 
two-step potentials. As for the HNC integral equation 
theory, it presents the best agreement in the region 1 < 
r/cr < Ai in the cases A-C, i.e., when ei > and £2 > 0. 
Even in those cases, however, the RFA is as accurate 
as or more accurate than the HNC theory in the region 
Ai < r/a < A2. For larger distances the RFA and HNC 
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FIG. 9. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system corresponding 
to case G (Ai — 1.25, A2 — 1.5, ei = 0, £2 = — e) at p* = 0.5 
and T* — 1.26193 with simulation results. Solid line: RFA; 
dashed line: HNC; dotted-dashed line: PY; symbols: REMC 
data. The inset shows the interaction potential and the tem- 
perature value. 
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FIG. 10. (Color online). Comparison of the different theo- 
retical approaches to compute the rdf of the system corre- 
sponding to case H (Ai = 1.25, A2 = 1.5, ei — e, €2 = —e) 
at p* = 0.5 and T* = 1.26193 with simulation results. Solid 
line: RFA; dashed line: HNC; dotted-dashed line: PY; sym- 
bols: REMC data. The inset shows the interaction potential 
and the temperature value. 
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FIG. 11. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system corresponding 
to case D' (Ai = 1.5, A2 = 2, ei = — e, £2 = e/5) at T* = 2 and 
(a) p* = 0.2, (b) p* — 0.4, and (c) p* — 0.75 with simulation 
results. Solid line: RFA; dashed line: HNC; dotted-dashed 
line: PY; symbols: MC data.— The inset shows the interac- 
tion potential and the temperature value. 



predictions are almost indistinguishable. In the rest of 
the cases (D-H), the PY and HNC curves are generally 
very similar, the best global performance being obtained 
with the RFA. 

As a final illustration, we present in Figs. [11] and 
1121 a comparison of the results we have obtained with 
the different theoretical approaches with those of MC 



simulations^ of two systems with an intermolecular po- 
tential of the type of case D (square well -I- square bar- 
rier), but this time having the following values of the 
parameters: Ai = 1.5, A2 = 2, ei = — e, £2 = e/5 (system 
D') and Ai = 1.5, A2 = 2, ei = — e, £2 = 2£/5 (system 
D" ), respectively. For those systems we have considered a 
fixed temperature T* = 2 and the three reduced densities 



"So 



bo 







FIG. 12. (Color online). Comparison of the different theoreti- 
cal approaches to compute the rdf of the system corresponding 
to case D" (Ai = 1.5, A2 = 2, ei = -e, t2 = 2e/5) at T* = 2 
and (a) p* = 0.2, (b) p* = 0.4, and (c) p* = 0.75 with sim- 
ulation results. Solid line; RFA; dashed line: HNC; dotted- 
dashed line: PY; symbols: MC data.— The inset shows the 
interaction potential and the temperature value. 



p* = 0.2, p* = 0.4, and p* ^ Q.lhM At the lowest density 
[p* — 0.2) the three theoretical approaches provide excel- 
lent results, with a slight superiority of PY and HNC. At 
the intermediate density {p* = 0.4), however, the RFA 
beats the PY and HNC results (which are practically in- 
distinguishable from each other). Finally, at the highest 
density [p* — 0.75) the RFA becomes clearly worse than 



the PY and HNC predictions (which are again hardly dis- 
tinguishable) , especially in the case of the higher barrier 
(caseD"). 



IV. CONCLUDING REMARKS 

In this paper we have presented a rather systematic 
study of the rdf of fluids whose molecules interact via a 
potential with a hard core plus two piece-wise constant 
sections of the same width and different heights (cither 
wells or shoulders), which include the RFA approach, our 
own REMC numerical experiments, and the numerical 
solution of the PY and HNC integral equations. We 
have considered eight representative classes of systems 
(see Figs. [1] and [2]) which cover all the possible topolo- 
gies of hard core plus two-step potentials. They include 
a purely repulsive potential (system A), potentials with 
an outer repulsive barrier and an inner attractive well 
(systems B-D), a purely attractive part outside the hard 
core (system E), and potentials with an outer attractive 
well and an inner repulsive barrier (systems F-H). Four 
of these systems (A and F-H) belong to the class of core- 
softened potentials. 

As Figs. IMTUl show, it is fair to state that for the fixed 
number density p* = 0.5 the agreement between the re- 
sults of our RFA formulation and those of the REMC 
experiments is very satisfactory in all instances. This is 
specially rewarding in view of the fact that the reduced 
temperature we have chosen to illustrate our findings 
(T* = 1.26193) is rather low and represents a stringent 
test of our theory. 

We have complemented our study with the comparison 
of the results we get and the MC data of Ref. [ij for two 
cases (labeled here as D' and D") in which the potential 
is of the form of the one of case D but with a wider width, 
a lower outer barrier, and a fixed reduced temperature 
T* = 2. Again, as observed from Figs. [TT] and [T^l the 
performance of the RFA approach is rather satisfactory, 
except at the highest density p* — 0.75. The tendency 
of the RFA to fail at high densities and wide potentials 
widths was already documented in the case of the pure 
square- well potential (with e2 = 0)i^ 

Concerning the comparison between our present ap- 
proach and the usual integral equation approach in the 
theory of liquids, we have seen that the RFA is rather 
simple, requires much less numerical labor (in these cases 
only the solution of two coupled transcendental equa- 
tions), captures correctly all the oscillations of the rdf, 
and is reasonably accurate. In fact, it is always superior 
to the PY equation [except for the cases of Figs. [TlT c) 
andfTWc)] and in most of the instances it is of comparable 
accuracy or better than the HNC equation. Hence, this 
constitutes further evidence of the usefulness of the RFA 
methodology that we have used for the computation of 
the structural properties of different hard-core fluids^S^ 

As a final point, we want to stress that the results 
of this paper, together with the earlier ones^ encour- 



age us to consider the problem in which the number of 
steps in the potential is much greater, leading in the limit 
ri — > cxD to the very interesting case of the structure in a 
fluid whose molecules interact via, for instance, a Jagla 
potential. We plan to undertake such a task in the future. 
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